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ABSTRACT 

We explore the observational characteristics of jet-driven supernovae by simulating bipolar-jet-driven ex- 
plosions in a red supergiant progenitor We present results of four models in which we hold the injected 
kinetic energy at a constant 10^' ergs across all jet models but vary the specific characteristics of the jets to 
explore the influence of the nature of jets on the structure of the supernova ejecta. We evolve the explosions 
past shock-breakout and into quasi-homologous expansion of the supernova envelope into a red supergiant 
wind. The simulations have sufficient numerical resolution to study the stability of the flow. Our simulations 
show the development of fluid instabilities that produce pristine helium clumps in the hydrogen envelope. The 
oppositely-directed, nickel-rich jets give a large-scale asymmetry that may account for the non-spherical exci- 
tation and substructure of spectral lines such as Ha and He I 10830A. Jets with a large fraction of kinetic to 
thermal energy punch through the progenitor envelope and give rise to explosions that would be observed to 
be asymmetric from the earliest epochs, inconsistent with spectropolarimetric measurements of Type II super- 
novae. Jets with higher thermal energy fractions result in explosions that are roughly spherical at large radii 
but are significantly elongated at smaller radii, deep inside the ejecta, in agreement with the polarimetric ob- 
servations. We present shock breakout light curves that indicate that strongly aspherical shock breakouts are 
incompatible with recent GALEX observations of shock breakout from red supergiant stars. Comparison with 
observations indicates that jets must deposit their kinetic energy efficiently throughout the ejecta while in the 
hydrogen envelope. Thermal energy-dominated jets satisfy this criterion and yield many of the observational 
characteristics of Type II supernovae. 

Subject headings: supernovae: general - hydrodynamics - instabilities - shock waves 



1. INTRODUCTION 

The observations of core-collapse supernovae are replete 
with evidence of asphericity. Hints of aspherical structure 
have been found in many supernovae (SNe) of various types, 
including both Type Ib/c (Maund et al.|2007a b ; Modjaz et al. 



2008 , Wang & Wheeler 2008 ) and Type II events (Fassiaetal. 



1998„Leonard et al..2001,.Elmhamdi et al.|2003„Chugai et al. 



2005 Leonard et al |2006 i. Challenges to the classic spheri 
cal model of supernovae (SNe) have been building since the 
observation of Type II SN 1987 A, wherein asphericity was 
inferred from net polarization measurements (Jeffery 199j]l 
and from subst ructure in spectral lines ( i .e., the 'Bochum' 
event, Hanuschi keraL][T95^ [Lucyl[T988| [PhilUps & Heafl>] 
cote 1989| ). Mixing one or two nickel clumps outward in 
the ejecta has been shown to reproduce the substructm^e of 
flie Hq line in SN 1987A dSpyromiUo eraL][T990) |Utrobin| 
eraL][T995] l. With the Hubble Space Telescope (HST), we 
have gained late-time, resolved images of the structure of the 
nascent SN r emnant that drama tically reveal aspherical, elon- 
gated ejecta (|Wang et al.|2 002[ l. 

Supernova 1987A, while especially observable, was a sin- 
gle event, and with an atypical blue supergiant progenitor; 
however, typical Type II-P supernovae (SNell-P), thought to 
arise from a red supergiant progenitor, show similar signs of 
asphericity and asymmetry. |Fassia et al.|(|T998) present in- 
frared spectra of Type II-P SN 1995V and discuss implica- 
tions of nonspherical elemental mixing. They find that in or- 
der to fit the He I 10830- A emission, significant amounts of 
^^Ni must be "dredged-up" to high velocities in order to ex- 
cite the He I line through radioactive decay. Fassia et al. also 
conclude that about 10% of the He mass must be contained in 
pure He clumps that are not microscopically mixed with hy- 



drogen. When mixed with hydrogen, collisional ionization of 
hydrogen (Penning ionization. B ell 1970, Chugai 19 9 1 ) is the 
dominant de-excitation mechanism for He in the hydrogen en- 
velope. This analysis suggests the need for a mechanism that 
produces significant outward propagation of nickel and the 
presence of fluid instabilities that lead to helium clumping. 

The detailed shape of spectral lines provides indirect evi- 
dence of asymmetry and asphericity of supernova ejecta. We 
can gain a direct measure of the shape of the supernova pho- 
tosphere and lin e-forming regions from continuum and line 
polarization (see Wang & Wheeler|2008 ). While the sample 
size of Type II-P SNe for which polarization measurements 
have been made is modest, the data to-date suggest asphericity 
which depends on the depth of the photosphere. Leonard et al. 
(2001) made spectropolarimetric observations of SN 1999em 
and find a net continuum polarization of about 0.2% on day 
7 after maximum light, increasing to 0.5% at day 165. Cit- 
ing the calculations of polarization from electron scattering 



spheroids by Hoflich (1991 1, Leonard et al. ( |2001) estimate 



the asphericity of SN 1999em to be ~7% on day 7 and ~10% 
on day 165. As they discuss, the late time polarization may be 
reduced due to the decreased optical depth to electron scatter- 
ing in the SN atmosphere, implying the intrinsic asphericity 
of SN 1999em may be greater than their estimates based on 
the electron scattering calculations of|Hbflich ( 1991 ). Further 
hints at the asphericity of the ejecta in SN 1999em are dis- 
cussed by Elmhamdi et al.| ( ["2003[ ). They report substructure 
in the Ha and He I lines that is reminiscent of the 'Bochum' 



event i n SN 1987 A. The spectral modeling of Elmhamdi et al. 
( 2003 1 shows that an asymmetric ^^Ni distribution may ac- 
count for the structure of these lines. 

A striking, and provocative, example of time-dependent 
SNell-P polarization is found in SN 2004dj. In this event, the 
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continuum polarization was practically zero, within the errors, 
during the plateau phase ( [Leonard et al.|2006) l, but during the 
transition from photospheric to nebular phases, as the photo- 
sphere receded to reveal the core of the SN, the polarization 
increased dramatically to 0.56% on day 91, implying an as- 
phericity of ~30%, and then slowly decreased until no polar- 
ization was detected on day 271. Leonard et al. ( 2006| l show 
that the decay of the polarization signal is roughly consistent 
with the expected decrease in electron scattering opacity of 
the receding photosphere. This seems to imply that the core of 
the explosion was significantly aspherical but was masked by 
a roug hly spherical envelop e. As discussed by [Leonard et al.| 
(2006 1, and first argued by Wang et al. (2001 1, when com- 
bmed with data from core-collapse SNe of other types (i.e., 
Ib/c), this may indicate that the core-collapse SN mechanism 
is similar across all spectral types and is intrinsically aspher- 
ical. Similar to SN 1999em, the spectral Hnes of SN 2004dj 
imply an asymmetric structure in the ejecta, inferred from the 
secondary peaks to the red and blue of the central peak. Spec- 
tral modeling including the energy deposition from radioac- 
tive decay of ^^Ni show that a bipolar nickel distribution in 
a spherical hydrogen envelope can reproduce the substructure 
of the Ha, [O I], and [Ca II] lin es due to asymmetri c line ex- 
citation from the decay of ^^Ni ( Chugai et al.|2005 1. Further- 
more, Chugai' ("2006 ) shows that a bipolar nickel distribution 
in a spherical envelope can also approximately account for 
the time evolution of the decreasing continuum polarization 
during the nebular phase. 

Based on the preceding discussion, explosion models for 
SNell-P must account for "dredge-up" of nickel to higher 
expansion velocities than expected from spherically sym- 
metric explosions, clumps of pristine helium in the hydro- 
gen envelope, an asymmetric and possibly bipolar distribu- 
tion of nickel, and roughly spherical electron-scattering pho- 
tospheres during the plateau phase that become remarkably 
aspherical as the photosphere recedes below the hydrogen 
envelope. The jet-driven model of core-collapse SNe (e.g.. 



Khokhlov et al.|19991[H5flich et al.|200Tl[Maeda & Nomoto[ 



2003| l has features that may produce some of the necessary 
observed characteristics. In this model, the explosion is driven 
by nonrelativistic bipolar jets that send strong bow shocks into 
the progenitor envelope. These bows shocks intersect near the 
equatorial plane of the star and establish an equatorial outflow 
that complements the polar outflows and completes the explo- 
sion. The material that comprises the jets is elongated along 
the jet axis and is expected to be nickel-rich owing to the jets' 
origin in the inner core of the explosion. Fluid instabilities at 
the interface between the helium core and hydrogen envelope 
may produce clumps of helium in the hydrogen envelope, as 
seen in non-jet-driven supernova calculations (e.g.,'Ki fonidis[ 
|H al. 2003 2006). The character of such an unstable flow, 
however, has not been investigated for the case of jet-driven 
SNe. It is also unclear that a jet-driven explosion in a red su- 
pergiant progenitor can produce a photosphere that is nearly 



spherical early-on and becomes aspherical later (see Hoflich 
l^al. 2001). 

In this work, we present a numerical study of the struc- 
ture of jet-driven SNe with adequate resolution to explore the 
nature of any unstable flows and time evolution sufficient to 
make a direct comparison with observations. We assume that 
bipolar jets are formed as a result of the core-collapse pro- 
cess and do not self-consistantly simulate the jet production 
process. Collapsing stellar cores could produce bipolar jets 
in the presence of rotation and magnetic field amplification 



( [Wheeler et aT:[[2000l |2002l [Akiyam a et al.|| 2003| l. The pro- 
cess of magnetic jet formation in collapsing stars has been 
studied in previous simulations ("LeBlanc & Wilson' 1970 Shi- 
bata et al. 2006; Obergaulinger et al. 2006a b; Burrows et al. 
|2007a( Komissarov & Barkov 2007), however, the rotation 
rate and initial magnetic field structure of the progenitor stars, 
two very critical properties, are not well understood from first 
principles. Therefore, we vary the parameters of the jets in 
four, physically-plausible jet models in order to explore the 
importance of the nature of the jets on the final structure of 
the SN. The progenitor model in all cases is a red supergiant 
star of 15 Mq, roughly analogous to Betelgeuse, and the mod- 
els are evolved to nearly six days after the initiation of the jets, 
several days after shock breakout. 

Many previous studies have focused on highly energetic jets 
in connection with Gamma-Ray Bursts (GRBs) (e.g., [Mac-| 
Fadyen et al.|2001 Zhang et al . 2003 ) a nd explosions o f Pop 



III stars (e.g., 
focus explicit 



Maeda & Nomoto,2003HTominaga|2007 |). We 
y on explosions of typical supernova energies 



(10^' erg) driven by nonrelativistic jets. We use a progenitor 
star of solar metallicity that is not extremely massive. The 
goal of this work is to explore the observational features of 
jet-driven supernovae in common Type II events. 

This work is organized as follows. In ^we describe our 
numerical approach and discuss the details of our jet models. 
In ^we provide a narrative of the general evolution of the ex- 
plosions, from shock propagation in the core to eruption from 
the progenitor's surface. The transport of nickel to radii and 
velocities greater than those of lighter elements is a feature of 
a jet-driven explosion model in which the jets contain a sig- 
nificant amount of nickel. In ^we discuss the significance 
of this source of overturn, as well as the development of fluid 
instabilities that lead to the mixing of elements, such as nickel 
and helium. In ^we discuss the resultant shapes of the ex- 
plosions and discuss the evolution of the surfaces of constant 
electron scattering optical depth. In ^we calculate the bolo- 
metric shock breakout light curves of our models and com- 
pare them to GALEX observations of supernova shock break- 
out from red supergiant stars. We discuss the implications of 
our study and list our conclusions in ^ 

2. THE SIMULATIONS 

We use the adaptive-mesh refinement (AMR) hydrodynam- 
ics code, FLASH ( Fryxell et al.[[2000| in 2D spherical -polar 
coordinates. FLASH is an Eulerian, higher-order Gudonov- 
type code that solves the hydrodynamic equations using a 
modified, piecewise-parabolic method (PPM) based on that 
of |Colella & Woodward| ( |1984) . We simulate a fluid with 8 
atomic species: H, He, C, O, Ca, Mg, Si, and Ni. Since in 
the present work we do not seek to precisely calculate result- 
ing metal abundances, we do not simulate the transmutation 
of species via nuclear burning. We use FLASH'S Helmholtz 
equation of state (EoS) that includes contributions to the in- 
ternal energy from ions (via the ideal gas law), electrons, 
positrons and radiation in thermal equilibrium. We calcu- 
late the gravitational potential from an evolving, spherically- 
averaged density profile. The refinement criterion was calcu- 
lated based upon the se cond spatial deriva tives of the density 
and pressure (see §2 of Fryxell et al.|2000 1 



We have run four different models with varied parameters 
for the bipolar jets. All models use a 15 Mq red supergiant 
progenitor (model "sl5s7b2" of Woosley & Weave r) 1995| for 
initial conditions. We map this model into 2D spherical geom- 
etry (r, 6) and place the initial inner boundary just outside the 
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iron core 3.8 x 10*^ cm, or ^ 1.6 Mq in mass); we do not 
simulate the collapse of the iron core. For all of our models, 
the angular extent of our 2D grid is a full 180°, correspond- 
ing to Ait steradians in 3D. We use seven levels of refinement. 
The base-level number of angular zones is 16 throughout the 
simulation, whereas the base-level number of radial zones is 
changed as the simulation proceeds to optimize the resolution 
(see re-gridding discussion below). The maximum angular 
resolution is 0.17°throughout the entirety of the simulations. 
All spatial derivatives of physical variables vanish at the in- 
ner and outer radial boundaries. Material is allowed to flow 
out of the computational domain through the inner and outer 
boundaries, but no material, other than that of the jets, is al- 
lowed to enter the domain through the inner boundary (the 
so-called "diode" boundary condition of Zin gale et al.|2002l l. 
Mass that flows out of the domain through the inner boundary 
is added to the central point mass contribution to the gravi- 
tational potential. The initial model is spherically symmetric 
and all velocities are initially set to zero. In the region beyond 
the surface of the progenitor (~ 3 x lO'-' cm), we place a sim- 
ple 1 /r^ density profile with a constant outflow velocity of 10 
km s"' that corresponds to a red sup ergiant wind with a mass- 
loss rate of lO"'* M© yr"' (see, e.g., 'Dwarka das|2005 1. Each 
model requires approximately 10,000 CPU hours to complete. 

In all cases, the north and south jets are identical in injec- 
tion method and initial parameters. We keep the total injected 
energy (10^' erg) and total time-integrated power (5 x 10^° 
erg s"') constant in all four models. We set the jets to be 
composed entirely of nickel. In reality, the composition of 
the jets may be time-dependent and not necessarily domi- 
nated by nickel. Simulating the jet formation mechanism self- 
consistently would be required to determine the exact compo- 
sition of the jets. The jets are injected into the progenitor 
model along both poles with an opening half-angle of about 
15° for a total of two seconds, roughly the free-fall time of the 
inner core. During the first half-second, the injection veloc- 
ity is kept fixed. After that, the velocity is linearly reduced to 
zero after two seconds. The time-dependent power of all jet 
models may then be described by 



f <0.5 s, 



P{t) = 



|(2-f)Po, 0.5s<f<2s, 



(1) 



f > 2 s, 



where is the initial jet power. The total energy (both ki- 
netic and internal) deposited by the jets in all models is 10^' 
erg, split equally between the two jets in each model. By set- 
ting the energy and power equal in all models, the jets are 
described by just two parameters: maximum injection veloc- 
ity and total injected mass. This is the origin of our naming 
scheming for the models: vAmS, where A is the velocity in 
units of 10^ cm s"' and B is the mass injected in units of hun- 
dredths of solar masses. 

Our method of jet injection is similar to that of |Khokhlov| 
[eTaL] ( [T999] l and [Hoflich et aT] ( [2001) except that by fixing 
the power of the jets as well as the total energy, the physical 
parameters are further constrained. We do not connect the 
accr etion rate into the inner boundary to the jet power (a s in, 
e.g.. 



MacFadyen et al. 2001 



nature of this accretion and t 



Maeda & Nomoto 2003 1, the 



le material that is accreted is 



lost in our calculations other than for its contribution to the 
gravitational potential. Notably, most of the silicon that lies 
just inside our computational domain falls through the inner 
boundary and so the final ejecta in our calculations have very 



little silicon. Since we are not conducting a careful study of 
nucleosynthesis and abundance yields this loss is acceptable. 

The parameters of the jets are given in Table [T] along with 
ambient values from the progenitor for comparison. Model 
v3ml2 is v ery similar to the jets employed by |Khokhlov] 
et al. ( 1999] ), however, in that work the progenitor was a he- 



lium star. The speed and density of v5m06 correspond ap 
proximately to those of the magnetically formed jets of |Bur- 
|rows et aL] ( |2007a| l. Models vlml2 and v6m04 are simi- 



lar in injection velocity to models simulated in Hoflich et al. 



(2001 1. While models v3ml2, v5m06, and v6m04 differ sig- 
nificantly in velocity and mass, the total energy in all three 
models is dominated by kinetic energy; hence, we shall re- 
fer to these models as the 'kinetic energy' models. This sets 
model vim 12 apart as the only model in which the thermal 
energy in the jet is greater than the kinetic energy. The ther- 
mal to kinetic energy fractions for the four models, found by 
integrating equation ([T]l, are 8 (vlml2), 0.05 (v3ml2), 0.06 
(v5m06), and 0.1 (v6m04). 

As the jets and shock waves expand to larger radii, we use 
a data regridding method to expand the grid and optimize the 
resolution, allowing us to achieve high-resolution, long time- 
scale simulations in reasonable amounts of computer time. 
While expanding the grid we also increase the radius of the 
inner radial boundary, excluding the small zones in the center 
of the grid where the Courant condition is very limiting and, 
thus, increasing the minimum time step. At each re-gridding 
we decrease the radial resolution but keep the angular reso- 
lution fixed throughout the entire simulation. We change the 
grid resolution parameters a total of ten times throughout the 
course of the simulations. The details of our changes in grid 
parameters are given in Table [2] All four models are evolved 
for a total 5.79 days, covering a spatial range from ^ 4 x 10^ 
cm to ~ 10'^ cm. 

3. DYNAMICS AND EVOLUTION OF JET-DRIVEN 
EXPLOSIONS 

In this section we provide a general analysis and discussion 
of the dynamics of our explosion models. The dynamics of 
the various jet models are primarily determined by the nature 
of the jets themselves, however, the structure of the progenitor 
star also affects the behavior of the blast waves. We show the 
radial chemical and density structures of the progenitor in Fig- 
ure[T] While the models we present here are multidimensional 
in nature and not well described by one-dimensional models, 
we can elucidate the dynamics by considering the evolution 
in the polar and equatorial regions separately. In Figure |2] we 
show the shock velocities in the polar and equatorial direc- 
tions as functions of radius in the progenitor star for models 
vlml2 and v3ml2. 

The evolution of the shock velocities in the polar and equa- 
torial directions is described qualitatively by the similarity so- 
lution for spherical adiabatic blast waves (Taylor 1946, Sedov] 
|1959[ [Chevalier 1976| l in power-law atmospheres described 
by p{r) (X r""^. Such solutions give the shock radius in time 
as Rshit) = /:i/(5-")f 2/(5-w)^ where is a constant, accleration- 
like term dependent on the central density of the atmosphere 
and the energy of the explosion. We can then write the shock 
velocity as a function of radius as 



Vsh(-'?sh) = Vsh.O 



Rsh 



sh.O 



cj/2-3/2 



(2) 



where Vsh.o is the shock velocity at some initial radius Rshfl- 
Equation ([2|i shows that for uj <3 (oj > 3) the shock deceler- 
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ates (accelerates), and for w = 3 the shock velocity is constant. 
In Figure[T|we also plot the power-law slope of the progenitor 
density profile, ui = -d\np/d\nR. Comparison with Figure l2l 
shows that, as expected, when > 3 the shocks accelerate and 
when uj <3, they decelerate. The polar shocks in all models 
depart from this behavior in the hydrogen envelope due to the 
growth of instabilities in the post-shock flow. We discuss how 
instabilities affect the shock evolution in more detail in ^ 

The remainder of this section is a basic description of the 
dynamics of our explosion simulations. It is included for 
completeness, however, the evolution is non-trivial and we 
attempt to desc ribe only what is most pertinent to the ensuing 



analysis. In 5 3.1 



we describe the evolution of the kinetic 

pro- 



energy-dominated jets in the progenitor core (j 3.1.1 
genitor envelope ( ^3.1.2| i, and circumstellar wind ( |3.1.3 1. In 



p.2| we give an analogous discussion for the thermal energy- 
dominated jets. In 5 3.3 we discuss the conditions that lead to 



instability growth that reaches the forward shock. Readers not 
interested in the detailed kinematics may skip to ^where we 
describe the development of material clumps in the jet mod- 
els. For reference, the transition from the C/O core to the He 
core is at a radius of about 2x10^ cm and the transition from 
the He core to the H envelope is located at about 6 x 10'" cm 
(see Fig. [T]). 

3.1. Kinetic Energy Jets 

The jet models experience three major phases of evolution: 
expansion through the core, expansion through the hydrogen 
envelope, and shock breakout and expansion into the circum- 
stellar medium. We illustrate these phases of evolution in Fig- 
ures|3]-|5]for model v3ml2. The two higher-speed kinetic en- 
ergy models evolve in a manner very similar to model v3ml2, 
and we will, therefore, restrict our detailed discussion in the 
following sections to v3ml2 alone. We give a comparison of 
the final density structures of these three models in Figure [6] 
The only significant difference in evolution as compared to 
model v3ml2 comes early while the jets are still in the pro- 
genitor core. In models v5m06 and v6m04, Kelvin-Helmholtz 
(KH) roll-ups grow much faster along the jet edges than in 
model v3ml2. This is attributed to greater shear with the 
progenitor structure in these models and, hence, a faster KH 
growth rate ( Chandrasekhar| 1 96T I . Jet model v3ml2 has pa- 
rameters taken directly from those of |Khokhlov et al.| ( |1999| l, 
except that our progenitor model has an intact hydrogen en- 
velope. We can compare our results with those of Khokhlov' 
|et al.| ( [l999j up to the point when the jets exit the helium core 
at around 6 x 10'" cm, which occurs at around 40 seconds, or 
the end of the first stage defined above. 

3.1.1. Evolution in the Core 

The injection of the jets into the core of the progenitor 
drives bow shocks that travel ahead of and laterally to the jets 
(see Fig. |3]l. Reverse shocks move away from the forward 
shocks. Contact discontinuities exist between the jet material 
and the core material swept up by the forward shocks. Vor- 
tices develop on the outer edges of the jet fluid (Figure [3] left 
panel). The edges of the jets themselves are also subject to 
small wavelength KH instabilities, but for the case of model 
v3ml2, these instabilities grow slowly and do not become dis- 
ruptive to the overall structure of the jets; this is in agreement 
with the findings of |Khokhlov et aL] ( |1999| l. For the other 



* The interested reader may view movies availabl e at jhttp : / / www . as ■ | 
[utexas ■ edu/$\sim$smc/SNmovies .htmll 



kinetic energy models, v5m06 and v6m04, the KH instabil- 
ities grow more quickly. After the jet injection stops at 2 s, 
the low-density region cleared out by the fast moving jet fluid 
closes, or heals, from the inner boundary out to the reverse 
shocks. This healing can be understood as the result of the 
low-pressure "vacuum" created by the cessation of the jet in- 
flow; the high pressure in the core of the star then acts to close 
this vacuum. Around this same time, about 10 - 15 seconds, 
the contact discontinuities show growth of Rayleigh-Taylor 
(RT) fingers (Figure |3] middle and right panels). Growth of 
these RT fingers is curtailed by shear with the surrounding 
post-shock flow. 

Along the equatorial plane, the two bow shocks cross one 
another, creating a dense, high-pressure "pancake" of material 
with net positive radial velocity that develops into an outward- 
moving torus of twice-shocked material (see Figure [3]l. This 
torus moves supersonically with respect to the surrounding 
progenitor and drives a shock wave out in the equatorial plane 
(see Fig. [3] right panel). Crossing the entire helium core takes 
approximately 40 seconds for the polar shocks and about 90 s 
for the equatorial shocks. 

3. 1 .2. Evolution in the Envelope 

The jets-torus structure that develops in the helium core is 
substantially altered by passage through the hydrogen enve- 
lope. Initially, the shocks accelerate down the steep density 
gradient between the helium and hydrogen (see Figs. [T]&[2]), 
but once in the < 3 hydrogen envelope, a reverse shock de- 
velops as the forward shock is decelerated. The reverse shock 
sweeps up the jet material and helium into a thin, dense shell. 
At the jet-star contact discontinuities, Richtmyer-Meshkov 
(RM) fingers develop, seeded by the pre-existing RT fingers 
of jet material (see middle pane of Figure|4|i. The RM fingers 
grow at roughly the free-fall speed and are able to reach the 
forward shock and perturb it (see ^ and 'Miles '2008). Be- 
fore being reverse shocked, the edges of the equatorial torus 
develop KH rollers (Figure]?] middle and right panels). The 
strong reverse shock, however, wipes out much of this detail 
as it sweeps up a thin shell. This shell is itself unstable and 
develops RT and RM fingers that become prominent as the 
explosion proceeds through the envelope. 

3.1.3. Evolution Following Shock Brealcout 

Figure|5]shows the evolution of model v3ml2 as the shocks 
exit the star and enter the circumstellar environment. Around 
3x10"* seconds after the start of the simulation, the polar 
shocks break out of the surface of the progenitor model (left 
panel of Figure [5]l. Due to slight differences between the evo- 
lution of the north and south jets, the northern shock breaks 
out slightly before the southern. The source of this north- 
south asymmetry is artificial. The build up of small numerical 
error after many time steps leads to non-uniform refinement in 
the the north and south. This, in turn, causes the north-south 
solutions to diverge in the subsequent evolution. The shocks 
accelerate down the steep density gradient in the outer en- 
velope and enter the lower-pressure circumstellar region (see 
Fig. [TJ. The shocks travel laterally and nearly encompass 
the entire perimeter of the star before the equatorial shock 
breaks out, roughly 10^ seconds after the polar shocks break 
out (middle and right panels of Figure |5]). The reverse shocks 
sweep up another thin shell of material that is RT unstable. 
The simulation ends at 5 x 10^ seconds (5.79 days), at which 
time the explosion is in approximate free expansion. The final 
kinetic energy of the ejecta is 5 x 10^" erg for model v3ml2. 
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whereas for models v5m06 and v6m04 the final ejecta kinetic 
energies are both ^ 7 x 10^*^ erg. 

3.2. Thermal Energy Jets 

Model vlml2 is the only model in which the initial energy 
budget of the jets is dominated by internal energy, as opposed 
to kinetic energy. Since we have fixed the energy and power 
to be the same in all four models, the slower speed of vlml2 
requires higher temperature to attain equivalent total power. 
Because not all of the internal energy in the jets is converted 
to kinetic energy, vlml2 produces a slightly under-energetic 
explosion; the ejecta posses total kinetic energy of ^ 3 x 10^° 
erg. 

3.2. 1. Evolution in the Core 

Following the relatively low-speed jet inflow, the thermal 
energy-dominated jet fluid quickly expands in all directions as 
the initial internal energy is converted to kinetic energy. The 
jet fluid attains a maximum radial velocity of about 3x10^ cm 
s~^ (see Figure[7]for the early velocity field in this model). The 
jets widen much more quickly compared to the other jet mod- 
els due to a higher degree of lateral adiabatic expansion (see 
left panel of Figure |8]l. After the jet inflow shuts off, around 
two seconds, the reverse shock moves quickly down the axis 
to the center of the star. The reverse shock closes the evacu- 
ated region opened by the jet at approximately four seconds. 
This is in contrast with the kinetic energy models wherein 
the jet cavity closes mostly from the bottom up. The less- 
directed nature of the energy deposition in this model leads 
to the spread of a significant fraction of the kinetic energy to 
large angles from the jets (see As a result, the equatorial 
outflow is larger and more energetic at early times than in the 
kinetic energy models. Around 50 s, the polar shocks exit the 
helium core, followed by the lateral shocks at ~80 s (middle 
and right panels of Figure |8]l. 

3.2.2. Evolution in the Envelope and Shock Breakout 

In the hydrogen envelope, the shock structure of model 
vlml2 becomes effectively spherical (see Figure |9]). A thin, 
dense shell develops behind the reverse shock that is, as in the 
kinetic energy models, RT unstable (Figure [9] middle panel). 
As the reverse shock hits the jet-star contact discontinuities, 
RM fingers develop but do not have enough time to catch the 
shock. As shown in Figures|9](right panel) and[TO] RT fingers 
form in the thin shell and grow to large sizes before the shock 
breaks out of the progenitor at ~ 10^ s. The shock departs 
the star roughly spherically and an unstable shell forms at the 
wind-induced reverse shock, analogous to the the kinetic en- 
ergy models (Figure |9j right panel). At the simulation's end, 
the model is approximately freely expanding. 

3.3. High-velocity Instabilities 

In the kinetic energy models, the shock is overtaken by 
downstream instabilities (see middle and right panels of Fig. 
|4]). This can occur for RT fingers in the fast-growth phase or 
RM fingers with Atwood number near unity (see Miles 2008 ). 
In model v3ml2 and the other kinetic energy models, RM fin- 
gers from the jet fluid contact discontinuity are indeed able to 
reach the forward shock in the hydrogen envelope. As the 
forward polar shock enters the H envelope, a reverse shock 
forms and travels downstream slowing the fast-moving post- 
shock material to a new post-forward-shock velocity. The 
contact discontinuity between the jet fluid and the helium has 



initial perturbations arising from earlier RT instability. Once 
reverse-shocked, the contact discontinuity becomes RM un- 
stable and the initial RT perturbations grow into RM fingers. 
In the frame of the contact discontinuity, the impulsive insta- 
bility velocity is ( |Youngs|1986| l 

vrm=^^— Af/ , (3) 

A P\ + P2 

where At/ is the change in velocity caused by the reverse 
shock, A is the wavelength of the initial perturbation, and aq is 
the amplitude of the initial perturbation. Prior to encountering 
the reverse shock, the contact discontinuity has perturbations 
characterized approximately by a^/X ^ 0.5. In the frame of 
the contact discontinuity, for model v3ml2 we estimate that 
the reverse shock causes an impulse |Af/| ~ 4 x 10^ cm s~'. 
The contact discontinuity has a post-reverse-shock Atwood 
number (/9i-/92)/(/0i+P2)"^ 3/4, andso vrm ^ 1 x 10^ cms"'. 
Transforming back into the lab frame, this implies that being 
reverse-shocked speeds up the initial perturbations in the con- 
tact discontinuity relative to its pre-reverse-shocked velocity 
by about 1 x 10 cm s"', meaning that the dense RM fingers 
are now moving at vo+ vrm ~ 1.2 x 10'' cm + 1 x 10^ cm ~ 
1.3 X 10^ cm s"' in the lab frame. This is approximately the 
speed of the forward shock at the time the reverse shock im- 
pacts the contact discontinuity. 

If drag on the RM fingers is negligible, they will move bal- 
listically with the speed given by equation As the shock 
encounters the shallow density profile of the hydrogen enve- 
lope (see Fig. [T]), it will slow below the speed of the RM fin- 
gers and they will overtake the shock. This is consistent with 
what is seen in our simulations. In Figure[2]the shock veloc- 
ity departs from the similarity solution given by equation Q 
at a radius of about 6 x 10" cm, where the RM fingers catch 
up. The new shock velocity is then set by the speed of the fin- 
gers which initially agrees with the estimate given above. The 
RM fingers and, hence, the forward shock do not continue at 
the constant free-fall velocity, but undergo deceleration due to 
drag and shock dynamics. 

Fast-moving RM fingers also develop at the jet-star con- 
tact discontinuity in model vlml2, however the separation 
between the contact discontinuity and the forward shock is 
much greater in this model (Fig. [5] left panel) and the fingers 
catch the shock at a much larger radius (^7x10'^ cm). Also, 
only one finger catches the shock in this model and it grows 
immediately along the axis. Growth of this finger may be 
artificially enhanced by the numerical coordinate singularity 
along the pole. 

4. CLUMPS IN JET-DRIVEN SUPERNOVAE 

4.1. Nickel Clumps 

Jet-driven supernovae have a robust mechanism that trans- 
ports heavier elements beyond lighter: the jets themselves. It 
is plausible that the jets contain a large fraction of nickel, or 
even, as we have assumed, are entirely nickel. This, then, 
permits the effective transport of nickel to velocities and radii 
where it can contribute to hydrogen and helium excitation ob- 
served at early times in SNe. In Figure 1 1 we show the final 



radial distribution of the masses of H, He, and Ni for the four 
models. As can be seen for model vlml2, a large mass of 
nickel has been carried out into the hydrogen envelope where 
it can excite the Ha and He I 10830A lines through radioac- 
tive decay. The presence of a bipolar arrangement of nickel 
clumps has been inferred from spectral modeling of SNell- 
P ( jElmhamdi et"iL][2003t |Chugai et al.|2005| |Chugai|[2006l ), 
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and our jet-driven models produce such a configuration (by 
Ansatz). We plot the Ni, H, and He mass distributions as func- 
tions of z-direction velocity (v^ = vv cos 9) in Figure 12 As can 



be seen, the nickel distribution of vlml2 consists of two op- 
positely directed clumps moving at a typical velocity of 1500 
km i"'. Each clump is composed of about 0.05 Mq, though 
this is certainly an upper limit because of our assumption that 
the jets are composed entirely of nickel. Figure [T3] shows a 
plot of the ^''Ni mass fraction for model vlml2 at the end of 
the simulation, illustrating the nickel clumps. 

The mass distributions for the kinetic energy models 
(v3ml2, v5m06, and v6m04) are quite different than that of 
model vlml2. The mass distributions for these models are 
given in Figures 11-12 The kinetic energy models drive 



a large fraction of the nickel mass beyond the hydrogen and 
helium. This may reduce the amount of nickel available to 
excite hydrogen and helium spectral lines. Only about half 
of the Ni mass of each model overlaps in radius with H and 
He, or about 0.05 Mq for vlml2, 0.03 Mq for v5m06, and 
0.02 Mq for v6m04. The existence of bare nickel beyond any 
hydrogen may also have implications for the gamma-ray sig- 
nal from such explosions, as a higher fraction of gamma-rays 
from the decay of nickel would escape. Figure [T2] shows that 
the bipolar nickel clumps reach much higher relative veloci- 
ties 5000 km s"') than for model vlml2. The final nickel 
distribution is shown for v3ml2 in Figure [T4| 

4.2. Helium Clumps 



Fassia et al. ( 1998 1 show that He I 10830A emission from 



SN 1995V requires the transport of nickel to high velocities 
and the presence of pristine helium in the hydrogen envelope. 
In our models, RT fingers that develop at the He/H interface 
produce pure-He clumps in the H envelope (see Figures[4]and 
|9j middle and right panels). In Figure [15] we show the He 
clump mass fraction xhc = A^He.ciumped/A^He for our four mod- 
els along with the H mass distribution. We consider a cell to 
be a pristine helium clump if it is more than 90% helium by 
mass. While the clumped fraction at certain velocities can be 
quite high (xHe.max ^ 0.1-0.2), the total fraction of He mass 
in clumps within the H envelope is only 1-2% across all four 
models. This is roughly a factor of ten less than what |Fas^ 



sia et al. (1998) find gives a good fit to the He I emission. 



We note, however, that our clumped fractions are lower limits 
due to numerical resolution. Any He clumps that fall below 
the minimum cell size become numerically mixed with other 
material. Despite the low clumped fractions, the strong He I 
emission may still be recovered from our models. We find a 
nickel mass at large radii that is a factor of two larger than that 
assumed by Fassia et al. ( |1998| l, and if this amount of nickel 
is included in the spectral modeling, the exposure of the He 
to a larger amount of Ni might enhance the He I emission to 
levels in accord with the observations. 

5. ASPHERICAL EXPLOSIONS AND IMPLICATIONS 
FOR SUPERNOVA POLARIZATION 

The current observational data show that SNell-P have 
time-dependent polarization signals, indicating that the 
shapes of SN photospheres change as the explosions expand. 
In this section we discuss the shapes of our simulated explo- 
sions and their implications for SNell-P observations. The 
final shapes of our models are controlled by the effectiveness 
of the jets to transport kinetic energy from the jet direction 
to more equatorial angles. In Figure 16 we show the angu- 



models. The kinetic energy distribution of model vlml2 at 
the end of the simulation is nearly flat, indicating a roughly 
spherically-symmetric energy distribution. The other three 
models leave significant amounts of kinetic energy in angles 
associated with the jets. This implies that the kinetic energy- 
dominated jets are "ballistic" and plow their way through the 
progenitor without losing much kinetic energy to more equa- 
torial angles through fragmentation and expansion. 

Carrying our simulations to quasi-homologous expansion 
allows us to explore the overall geometry of the explosions 
at times when supernovae are actually observed. As dramat- 
ically demonstrated by [Leonard et al.| ( |2006| l, there is com- 
pounding evidence that SNell-P exhibit very low polarization 
through the plateau phase but, as the photosphere recedes be- 
low the hydrogen envelope, the polarization increases greatly. 
Since the continuum polarization is primarily the result of 
electron scattering, we may compare the shapes of surfaces 
of constant electron scattering optical depth, Tes, in our mod- 
els with SN geometries inferred from polarization measure- 
ments. To compute Tes, we first expand our simulation results 
to the desired time using the relation for homologous expan- 
sion, r = vt. We calculate the new temperature of our models 
by assuming that the expansion is adiabatic. Considering only 
radiation and gas pressure, the internal energy is then 



flradT" 3 ksT 
€ = 1 = const. 

P 2 /iOTH 



(4) 



Assuming adiabatic expansion neglects radiative losses, how- 
ever, this will be adequate for qualitatively estimating the pho- 
tosphere shapes. The new gas temperature is found by solv- 
ing equation Q, which is then used to calculate the electron 
fraction arising from thermal ionization of hydrogen, xj, us- 
ing Saha's equation. We include contributions to the electron 
fraction due to gamma-ray ionizations from t he radioactive 



decay of ^^Ni in a manner similar to that of [Utrobin et al. 
( 1995 1. We first calculate the local gamma-ray luminosity on 
a cell-by-cell basis as 

Lm = (6.456-'/**-^ ''^y + 1.45e-'/'"-^ ''"'') x lO'^^mNi/M©), (5) 

where is the mass of nickel in a given cell. Assuming that 
the gamma-ray luminosity contributes only to hydrogen ion- 
ization, and entirely locally within a cell, we find the electron 
fraction due to nickel decay by iteratively solving 



a32xt^niidV = (1 -x"'^**)^— 

-^13 



(6) 



for x^, the electron fraction arising from the gamma-rays, 
where q;32 = 3.37 x lO"'-' cm-' s"' is the de-excitation rate at 
5000 K, rtu is the hydrogen concentration, dV is the cell vol- 
ume, 7/ = 0.39 is a hydrogen ionization branching ratio, and 
£13 = 12. 1 eV is the transition energy from the second excited 
state to the ground. Assuming that all gamma-ray energy is 
locally deposited within a cell is valid because, in our case, the 
typical gamma-ray mean free path, / = l/n^p, is small com- 
pared to the simulation cell size for a gray gamma-ray opacity 
of ^ 0.03 cm^ g~' due to elec tron scattering and absorp- 
tion ( Sutherland & Wheeler] 1984| . The total electron fraction 
is taken to be x = max[ l, xr+x^]. The electro n scattering 
optical depth is then (see Jeffery & Branch|1990| l 

f« 

:0.4 / xpudr, (7) 



lar distribution of kinetic energy at various times for the jet 



where is the hydrogen mass density and r is along the ob- 
server's line of sight. 
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In Figure 17 we show the evolution of the electron scatter- 
ing surfaces as a function of time for model vlml2, calculated 
using the method described above. Early-on, about 5 days af- 
ter shock breakout, the photosphere (res ~ 1) is slightly oblate. 
This is a result of the slightly aspherical breakout. The polar 
shocks erupt from the star about a third of a day before the 
equatorial shocks and have that much time to expand later- 
ally, and thus the kinetic energy of the polar shocks is diluted 
and the shock slows. On day 50, the photosphere is roughly 
spherical, in agreement with the early, plateau-phase polar- 
ization measurements of SN 2004dj (Leonard et al. 2006 ). As 
the photosphere recedes, its shape elongates along the jet axis. 
By day 100, the photosphere is a prolate spheroid of axis ra- 
tio ~ 1.25. This is about the epoch associated with the end 
of the plateau phase in SNell-P and about the time Leonard] 
|et al.| ( |2006 ) saw a dramatic rise in the polarization signal 
from SN 2004dj. Leonard et al. (2006i estimate an axis ratio 



of about 1.4 for SN 2004dj based upon the electron scattering 
spheroid polarization calculations of |Hoflich| ( |1991| l. Given 
the many uncertainties in the interpretations of both the ob- 
servations and our calculations, there is striking agreement 
between the photospheric evolution of model vlml2 and that 
of SN 2004dj. 

The kinetic energy models are far more asymmetric at all 
times, as is demonstrated for model v3ml2 in Figure 18 Axis 
ratios for this model range from 2 to 3 from the earliest epochs 
until the explosion begins the transition to the nebular phase. 
The shape of the photosphere on day 300 is actually more 
spherical compared with earlier times. The dramatic aspheric- 
ity of model v3ml2 is inconsistent with spectropolarimetry 
of SNell-P ([Leonard et al.||2001| [20061 [Wang & Wheeler 



2008| l. The shapes of the photospheres are similar for models 



v5m06 and v6m04. This seems to indicate that if SNell-P are 
driven by bipolar jets, such jets must be "stalled" by passage 
through a thick hydrogen envelope, meaning that they must 
efficiently lose kinetic energy to the envelope driving a more 
nearly spherical explosion. A similar conclusion was reached 
by Hoflich et al. ( 2001) ). Our thermal energy-dominated jet 
model satisfies this criterion. 

6. SHOCK BREAKOUT 

Recently, the GALEX satellite has observed two SNell-P 
at, or near, the moment of shock break out (jSchawinski et al.] 
[2008|[Gezari et al.|2008[ ). [Gezari et al.[ ( [2008| l model the hght 
curve as a spherical breakout, but as we show here, the break- 
out light curve will vary according to the asphericity of the 
shock breakout. To illustrate this, we have computed break- 
out light curves of our models by first finding the photosphere 
(Tes = 1; via the method outlined in ^ and then calculat- 
ing the bolometric blackbody emission at the photospheric 
temperature in a piecewise fashion across the photosphere 
and summing these blackbody contributions over the photo- 
sphere. The resultant light curves, along with the angled- 
averaged photospheric temperature, are shown in Figure [19] 
The significantly aspherical breakout of model v3ml2 is ap- 
parent in the "double-peaked" light curve. The breakout of 
the jets around 0.5 days, measured in time from the beginning 
of the simulation, is followed by a period of near-constant 
luminosity associated with the expansion of the jet cocoons. 
The equatorial shock breakout leads to the second peak at 1 .5 
days. The light curve of model vlml2 is not significantly 
different from that of a spherical shock breakout model (see 
Figs. 1 and 2 of Gezari et al.|2008f The slightly aspherical 
breakout of model vlml2 is evident in the light curve, but 



the first peak in this case is much smaller and is much closer 
in time to the large peak than in model v3ml2. Consider- 
ing the GALEX observations, the early peak and plateau in 
the light curve of model v3ml2 would have been detected for 
SNLS 04D2dc if the breakout were strongly aspherical. The 
absence of this early peak and plateau in the breakout light 
curve indicates the breakout was indeed nearly spherical. We 
do note, however, that the temporal resolution of the GALEX 
data is roughly one hour, such that the early peak in the light 
curve of v3ml2 would have been only marginally detected. 
The light curve of model vlml2 is in much better agreement 
with the GALEX data. While our method of computing the 
light curves is crude and neglects many important aspects of 
shock breakout (such as the radiative precursor), our theoreti- 
cal light curves do illustrate the significant differences that can 
arise from a strongly aspherical breakout and indicate that the 
characteristics of such a breakout could have been detected by 
GALEX, if they were present. 

7. DISCUSSION AND CONCLUSIONS 

We have performed 2D axisymmetric, high-resolution sim- 
ulations of stellar explosions driven by bipolar jets. We have 
followed the evolution of four different jet models to near free 
expansion. The energy budgets of three of our jet models 
(v3ml2, v5m06, and v6m04) are dominated by kinetic energy 
whereas that of one jet model (vlml2) contains a large frac- 
tion of thermal energy, and all jet models inject a total of 10^' 
erg. This distinction in the four models drastically affects the 
resulting structures of the models. The kinetic energy models 
remain highly elongated as they expand into the wind sur- 
rounding our progenitor model, whereas model vlml2 pro- 
duces an effectively spherical outer envelope with large de- 
partures from symmetry in the core. These shapes correspond 
to how efficiently the jet kinetic energies are isotropized. Only 
model vlml2 results in an effectively spherically-symmetric 
kinetic energy distribution (see Fig. 16 1. An anisotropic dis- 
tribution of kinetic energy may have implications for the mea- 
surement of supernova energies from observations. While ra- 
diative transport calculations of our results would be required 
to reproduce typical spectral lines used in estimating super- 
nova energies, it is likely that the kinetic energy models, with 
their asymmetric energy distributions, would yield viewing 
angle-dependent kinetic energy estimates. If a spherical ex- 
plosion model were assumed when analyzing the observations 
of an aspherical supernova, errors in the estimate of the kinetic 
energy could arise. 

We resolve the growth of fluid instabilities in the explo- 
sions. We show that the jet fluid is unstable and develops 
several RT fingers that, when reverse-shocked in the H enve- 
lope, serve as initial perturbations for rapid RM growth. The 
RM fingers travel at nearly the free-fall velocity and are able 
to overtake the shock in the kinetic energy models. The edges 
of the jets are also KH unstable, but the growth rate depends 
sensitively on the speed and density of the jet fluid. Model 
v3ml2 does not show significant KH growth, but all three 
other models do. The equatorial outflow in all four models is 
also subject to KH growth along the top and bottom surfaces, 
as predicted by [Wheeler et al.[([2008|l. 

In comparison with Khokh ov et al. ( 1999) 1, model v3ml2 
shows a similar large-scale structure up to the time the shocks 
depart the helium core (see their Figure 2). The extension in 
the jet direction is similar and the development of the equa- 
torial outflow is comparable. In contrast, the higher reso- 
lution employed in our model (Armin ~ 3 x 10^ cm versus 
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Ar^in ~ 4 X 10^ cm in Khokhlov et al.) allows us to see the 
development and growth of fluid instabilities, particularly in 
the jet fluid (see the right panel of Figure[3]l. 

We have compared the late-time geometries of our models 
with observations by calculating the electron scattering opti- 
cal depths under the assumption of adiabatic expansion (see 
Figs. [17] and 18 1. We show that the kinetic energy-dominated 
models (v3ml2, v5m06 and v6m04) are too aspherical at 
early times to match polarimetry of SNell-P. Model vlml2 
produces a photosphere that is nearly spherical for the first 
50 days of expansion and increases in prolate elongation to 
axis ratios t hat are in close agreem ent with those derived for 
SN 2004dj ( [Leonard et"aL] [2006) 1. Our simpUfied calcula- 
tion of optical depths neglect several important radiative ef- 
fects, however a full radiative transfer calculation is beyond 
the scope of this work and is unlikely to qualitatively change 
our results. 

The development of instabilities may be important to the 
overall understanding of core-collapse supernovae in a num- 
ber of ways. Substantial observational evidence shows that 
the transport of heavier material to radii ahead of lighter mate- 
rial is a ubiquitous phenomenon in core-collapse supernovae, 
especially SNell-P. This "overturn" of the original "onion- 
skin" configuration of the progenitor has been explained in 
neutrino-driven explosions by the growth of fluid instabilities 
(see |Kifonidis et al.|[2003| 2006) . The substructure of cer- 
tain strong lines (such as Ha) indicates large, fast-moving 
clumps of nickel in the ejecta. Similar to the calculations 
of |Kifonidis et al.] ( |2003[ ), many of the instabilities that de- 
velop in our models are of very small scale and it is unlikely 
that such small instabilities can explain the substructure of 
spectral lines, however, the jets are effectively a large, single- 
mode asymmetry capable of producing high-velocity bipolar 
clumps. The jet fluid, which we assume to be nickel, is trans- 
ported ahead of lighter materials, providing significant "over- 
turn." In our models, a large amount of nickel mass is trans- 
ported to radii where it can excite the Ha and He I lines, as is 
required by spectral modeling of SNell-P (jFassia et al.'1998S 
Elmhamdi et al.|2003[[Chugai et al.|2005t |Chugai 2006). Ki- 



fonidis et al.|p006| l show that the mode of the instabilities that 



grow in neutrino-driven explosions is sensitive to the nature of 
the neutrino emission. In jet-driven models, the jets create a 
robust, large-scale asymmetry that is not very sensitive to the 
detailed nature of the jets. 

We show that He clumps form due to RT instabilities in the 
transition layer between the H and He in the explosion ejecta 
in all four jet models. These clumps may contribute the pris- 
tine helium required by Fassia et al. ( 1998) to fit the strength 
of the He I 10830 A line, though the amount of clumping we 
find is below that required by Fassia et al. As we have noted, 
the extra nickel at high velocities in the jet models (about a 
factor of two greater than that assumed by Fassia et al.) may 
compensate for the low helium clump fraction by providing 
excess gamma-ray excitation. 

Evolving our calculations beyond shock breakout allows us 
to study the breakout characteristics of the explosion models. 
In Figure [T9[ we present the bolometric shock breakout light 
curves for models vlml2 and v3ml2. The nearly-spherical 
breakout of model vlml2 produces a light curve that is very 
similar to spherically-symmetric model light curves. The dra- 
matically aspherical breakouts of the kinetic models show a 
first peak in the light curve associated with the breakout of the 
jets followed by a roughly constant luminosity period created 



by the expanding and cooling jet "cocoons" and then a large 
peak as the equatorial shocks break out. The early peak and 
plateau in the light curves of the kinetic models would have 
been evident in GALEX observations of SNLS 04D2dc, if they 
were present. Their absence argues against the likelihood of a 
significantly aspherical break out for this SNell-P, as does the 
comparison to polarimetry. The light curve of model vlml2 
agrees with the observations for shock breakout for the two 
SNell-P ca ught by GALEX ( [Schawinski erai:][2008t |Gezari| 
let al.|2008|). 

In Table [3j we summarize the observational characteristics 
of SNell-P and the four jet models. We also include the char- 
acteristics a neutrino-driven SN model based loosely on that 
of Kifonidis et al. (2003 20061 for comparison. A multi- 



dimensional neutrino-driven model that includes fluid insta- 
bilities can be expected to account for transport of heavy el- 
ements to high ejecta velocities as well as the formation of 
pristine helium clumps in the hydrogen envelope; however, it 
is not clear that neutrino-driven models can consistently pro- 
duce fast-moving, bipolar nickel clumps in the ejecta. Like- 
wise, it is unlikely that these models can account for the time 
dependence of SNell-P polarization. Significantly aspherical 
neutrino-driven explosions may result from the standing ac- 
cretion shoc k instabihty (SASI; |BEndin et al. 2003 , Marek M 
'Janka'2007) or the acoustic mechanism (Burrows et al.|2006f 
,2007b) , however the large-scale structure of supernovae re- 
sulting from these mechanisms has not been explored. As 
we have discussed, the kinetic energy-dominated jet models 
produce explosions that are too aspherical at early times and 
are, thus, incompatible with polarization and shock break- 
out observations of SNell-P. In contrast, our thermal energy- 
dominated model (vlml2) shows good agreement with all of 
the observations we have discussed here. 

The comparison of observational evidence with our simu- 
lations consistently indicates that jets capable of retaining a 
large fraction of their initial kinetic energy while crossing the 
hydrogen envelope of a red supergiant star are unlikely to be 
present in SNell-P explosions. This is evidence that if jets are 
present, or are the driving mechanism in core-collapse super- 
novae, they must be jets that are capable of transferring their 
kinetic energy to the hydrogen envelope of a RSG. This may 
not require thermal energy-dominated jets, as we have pre- 
sented here, but it does seem to rule out fast, dense jets that 
remain ballistically detached from the surrounding ejecta and 
punch out of the progenitor with a large fraction of their initial 
kinetic energy intact. |H6flich et al.|(|2001|l and [Khokhlov & 



Hoflich (2001) have noted that if jets efficiently penetrate the 
core of a stripped envelope supernova, but are mostly trapped 
in the envelope of extended hydrogen-rich supernovae, then 
average stripped-envelope supernovae should have lower ki- 
netic energies than average SNell-P. Since this is not the case, 
Hoflich deduced that any jet energy must already be fairly effi- 
ciently deposited in the cores of supernova progenitors by dint 
of any jets being rather wide and slow. The conclusions of our 
current study, that thermal jets are favored over kinetic energy 
jets, is roughly in accord with Hoeflich's supposition. We 
note, however, that even our thermal jets that ultimately de- 
posit their energy rather uniformly in the hydrogen envelope 
are less efficient in depositing energy in the helium core. The 
energy deposition cannot be completely isotropized in the he- 
lium core without violating the observations from spectropo- 
larimetry that the inner core is substantially aspherical. How 
these two constraints of energetics versus asymmetry play off 
against one another needs to be studied in more detail. 
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TABLE 1 
Parameters of jet models 



Jet Model 


P (g cm 3) 


Tmax (K) 


P (erg cm ^) 


Vmax (cm/s) 


Mto, (Mq) 


v3ml2 


6.5 X 10^ 


2x 10' 


1023 


3.2 X 10' 


0.12 


v5m06 


2.2 X 10^ 


2x 10' 


1023 


5 X 10' 


0.06 


vlml2 


2x 10^ 


4.4 X 10' 


1024 


10' 


0.12 


v6m04 


105 


2x 10' 


3 X 1022 


6x 10' 


0.04 


ambient" 


8 X 10' 


2x 10' 


1.4 X 1023 







" The ambient values given are those of the progenitor model and the inner radius of 
our simulations: 3.8 X 10^ cm. 



TABLE 2 

Simulation stages and grid parameters 



Stage 


t, (s) 


tf (s) 


nn (cm) 


rout (cm) 




Ar„,i„ (cm) 










5 


3.82 X 10** 


3.2 X lO'O 


192 


2.6 X 10^ 


1 




5 


25 


3.82 X 10** 


7.5 X 10'" 


192 


6.1 X 10* 


2 




25 


100 


1.0 X 10' 


2.5 X 10" 


192 


2.0 X 10' 


3 




100 


300 


2.0 X 10' 


9.0 X 10" 


192 


7.3 X 10' 


4 




300 


1 X 103 


4.0 X 10' 


3.0 X 10'2 


320 


1.5 X 10* 


5 




1 X 103 


3x 103 


1.0 X 10'° 


6.0 X 10'2 


320 


2.9 X 10** 


6 




3x 103 


1 X lO"* 


2.0 X lO'O 


2.0 X 10'3 


352 


8.9 X 10** 


7 




1 X 10-* 


3x 10" 


5.0 X 10'" 


6.0 X 10'3 


352 


2.7 X 10' 


8 




3 X 10* 


1 X 10' 


l.Ox 10" 


2.0 X 10"' 


384 


8.1 X 10' 


9 




1 X 10' 


2x 10' 


2.0 X 10" 


4.0 X 10'" 


384 


1.6 X 10'" 


10: 


2x 10' 


5 X 10' 


5.0x 10" 


l.Ox 10" 


384 


4.1 X 10'" 



" The number of radial zones at the lowest level of refinement. 
^ The radial resolution at the highest level of refinement. 



TABLE 3 

Observational features of supernova models 



Observed characteristic of SNell-P 


vlm2 


v3ml2 


v5m06 


v6m04 


neutrino-driven" 


early polarization small 


yes 


no 


no 


no 


yes 


late polarization large 


yes 


yes 


yes 


yes 


no 


spectral substructure due to Ni clumps 


yes 


yes 


yes 


yes 


no'' 


pristine He clumps in H 


yes 


yes 


yes 


yes 


yes 


GALEX breakout light curve 


yes 


no 


no 


no 


yes 



" Based loosely on the results of |Kifonidis et al.H2003||2006^ . 



jKifonidis et al. (20061 do find large-scale nickel clumps that may account for spectral substruc- 
ture, however this result seems to depend sensitively on the nature of the neutrino emission. A 
better understanding of neutrino physics is needed in order to know if this result can be obtained 
consistently in all SNell. 
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E (cm) R (cm) 



Fig. 1. — Radial abundance (left) and density (right) profiles of the red supergiant progenitor model used as initial conditions for our simulations. In the left 
panel, X, is the mass fraction of the atomic species listed in the legend. Along with the density profile (sohd line) on the right, we plot the logarithmic slope of 
the density, ti; = —d\np/d\i\R (dashed line). The spikes in lu correspond to layer transitions in the progenitor: C/O-He (2 X lO' cm), He-H (6 X lO'" cm), and 
H-wind (3 X 10'' cm). 




10» 



10'" 



10" 10'^ 
R (cm) 



Fig. 2. — Shock velocity as a function of radius for models vlml2 (thermal energy-dominated, left) and v3ml2 (kinetic energy-dominated, right). The solid 
line is the shock velocity along the pole (6 = 0°) and the dashed lined is the shock velocity along the equator (d = 90°). For model vlml2, the shock structure 
becomes roughly spherical in the H envelope; the solid and dashed lines overlap until an instability finger grows up the axis and overtakes the polar shock around 
5 X 10'^ cm. The shock structure never becomes spherical in the case of model v3ml2 and instability fingers overtake the shock much earlier, around 7 X lO" 
cm. The equatorial shocks are well-defined only at radii greater than about 2 X 10' cm. 
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Fig. 3. — Snapshots of the density evolution in kinetic energy-dominated model v3ml2 before leaving the helium core of the star, which has a radius of 
~ 5 X lO'" cm. The axis lengths are given in centimeters and the density is plotted in g cm"^^. The simulation times are, from left to right, 2 s, 10.6 s, and 40 s. 
Note the stripping of jet material caused by shear with the stellar material and the subsequent development of small-scale, complex structure in the right panel. 
The three panels also show the development and evolution of the equatorial structure resulting from the collision of the jet-induced shocks at the midplane. The 
middle and right panels also show the beginning growth of RT fingers at the interface between the jet material and the shocked star material. 




R (xl0'12 cm) R (xl0~ 12 cm) R (xl0'12 cm) 



Fig. 4. — Snapshots of the density evolution in kinetic energy-dominated model v3ml2 while transiting the stellar envelope. The axis lengths are given in 
centimeters and the density is plotted in g cm"^. The simulation times are, from left to right, 700 s, 2200 s, and 26,000 s. In the equator, the "pancake"-like 
outflow has developed KH roU-ups as it has sheared through the envelope of the progenitor The leading KH roll-up detaches from the bulk of this flow and 
impinges on the equatorial shock, as shown in the middle panel. The RT fingers of jet fluid are swept up by the reverse shock into a thin, dense shell. The fingers 
serve as an effective initial perturbation to instability growth in this shell that, as seen in the middle and right panels, is able to reach up to and ahead of the 
forward shock, thus perturbing the forward shock structure. 
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Fig. 5. — Snapshots of the density evolution in kinetic energy-dominated model v3ml2 after shock breakout from the star. The axis lengths are given in 
centimeters and the density is plotted in g cm"^. The simulation times are, from left to right, 52,000 s, 130,000 s, and 500,000 s. This figure shows the significant 
lag between the breakout of the polar shocks and that of the equatorial shock. A thin, unstable shell of wind material is swept up by the forward shocks, analogous 
to that shown in Figure |4] The final structure (right panel) is significantly prolate, inconsistent with early polarization measurements of SNell-P. Note that the 
slight north-south asynmietries are artifacts of numerical inaccuracies and are amplified by shock breakout. 




Fig. 6. — The density structures of the three kinetic energy-dominated models v3ml2, v5m06, and v6m04 (left to right), at the end of the simulations, 5.79 
days. The final structures are remarkably similar given the differences in jet injection speeds and jet masses. 
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Fig. 7. — The radial velocity field of thermal energy-dominated model vlml2 at 1 s. The maximum velocity attained of ~ 3 X lO' cm s ' is much greater than 
the input velocity of lO' cm s"' due to the conversion of thermal energy to kinetic energy in the jet fluid. 
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Fig. 8. — Snapshots of the density field for thermal energy-dominated model vlml2 before leaving the helium core. The axis lengths are given in centimeters 
and the density is plotted in g cm"^^. The times of the snapshots are, from left to right, 2 s, 20 s, and 60 s. The left panel shows the rapid lateral expansion 
of the jet structure in this model. Shown also are large vortices that develop into turbulent structures, seen in the middle and right panels. The right panel 
shows the beginning growth of RT fingers in the interface between the jet and the shocked star material. The equatorial structure is apparent in the middle panel 
and is significantly broader than in the kinetic energy models. The edges of the equatorial structure are less sharp than in the other models, but still subject to 
short-wavelength KH instability (see Fig. |9}. 



COUCH, WHEELER, & MILOSAVLJEVIC 



15 




K (xl0~9 cm) R (xlO'lZ cm) R (xl0'12 cm) 



Fig. 9. — Snapshots of the density field for thermal energy-dominated model vlml2 in the progenitor envelope. The axis lengths are given in centimeters and 
the density is plotted in g cm"'. The times of the snapshots are, from left to right, 500 s, 6500 s, and 85,000 s. The left panel shows that the equatorial outflow is 
much broader than in the kinetic energy models and is still subject to the KH instability along its surfaces. As also seen in the kinetic energy models, the reverse 
shock sweeps up a thin, RT unstable shell that wipes away much of the complex structure that has developed up to this point. The right panel shows the RM 
finger of jet fluid that grows up the axis and catches the shock. 




J! (xl0'12 cm) R (xl0'12 cm) R (xl0'12 cm) 



Fig. 10. — Snapshots of the density field for thermal energy-dominated model vlml2 as the shock exits the star. The axis lengths are given in centimeters and 
the density is plotted in g cm"'. The times of the snapshots are, from left to right, 130,000 s, 155,000 s, and 500,000 s. This figure shows that the lag between 
polar and equatorial shock breakouts is much smaller than in the kinetic energy models. The resulting structure is nearly spherical. As in the case of model 
v3ml2, the reverse shock arising from shock breakout into the wind sweeps up a dense shell that is RT unstable. 
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Fig. 11. — Mass distributions versus radius for the four explosion models at the end of the calculations (5.79 days). The three kinetic energy models drive 
nickel ahead of the hydrogen and helium while the thermal model, vlml2 (upper left), ends with a nickel distribution very similar to that of hydrogen. Driving 
nickel beyond the hydrogen and helium may reduce the nickel available to excite spectral lines via radioactive decay. 
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Fig. 12. — Mass distributions versus z-velocity (Vj = VrCOsfl) for the four models. The oppositely-directed, bipolar clumps of nickel are evident. The kinetic 
energy models yield nickel clumps with typical velocities of 5000 km s~' whereas the clump velocities in model vlml2 are ~ 1500 km s~' . 
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Fig. 15. — Distribution of the clumped He mass fraction (XHe = 'WHc^clumpcd/^He) and H mass distribution in the four models. In all models, a small fraction 
of the He mass in the outer hydrogen envelope is clumped, [passia et al.|{l998j find that a much higher clumped fraction is needed to reproduce the He I 10830A 
line in SN 1995V; however, their modeling did not account tor sigmhcant nickel clumps at high velocities (see figs. |l 1| - |12| . 
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Fig. 16. — Angular distribution of kinetic energy for the four models at various times. Note that the thermal energy-dominated model (vlml2, upper left) has 
distributed its energy nearly uniformly in angle by the end of the simulation whereas the three kinetic energy-dominated models (other panels) retain significant 
kinetic energy near the jet axis. 
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Fig. 17. — Surfaces of constant electron scattering optical depth r^s using a constant opacity for thermal energy-dominated model vlml2 at several epochs: 6 
(end of simulation), 50. 100, and 200 days. The optical depths were calculated for a line-of-sight in the equatorial plane and the levels shown are Tcs = 1, 10,30. 
On day 6, the photospheric shape is slightly oblate with axis ratio of ~ 0.95, but by day 50 the photosphere is essentially spherical. On day 100, around the time 
of the observed jump in polarization in SN 2004dj, the axis ratio of vlml2 is ~ 1.25 and the photosphere is prolate. On day 200, the elongation of the inner core 
is apparent and the contour for Tcs = 30 is absent as the model is making the transition to the nebular phase. 




ii(10"cm) B(10'''cm) ii (10'" cm) H (10'' cm) 

Fig. 18. — Same as Fig. |17| except for kinetic energy-dominated model v3ml2. The shape of the photosphere is remarkably asymmetric at all times, in 
disagreement with polarimetry of SNell-P. 
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Fig. 19. — Bolometric shock breakout light curves for models vlml2 and v3ml2. The light curve of thermal model vlml2 is similar to that for a spherically- 
symmetric breakout and is consistent with GALiJX observations of shock breakout from SNLS 04D2dc. The kinetic model, v3ml2, has a double-peaked breakout 
light curve that would have been evident in the GALEX observations, but was absent. We also plot the angle-averaged photospheric temperature. 



